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THE STABILITY OE. LAMINAR ELOW PAST A SPHERE* 

By J . Pretsch 

summary' 



As a contribution to the problem of turbulence on a 
surface of rotation, the method of small oscillations is 
applied to the flow past a sphere. It was found that the 
method developed for two-dimensional flow is applicable 
without modifications. The frictional layer in the vicin- 
ity of the stagnation point of a surface of rotation is 
less stable against small two-dimensional disturbances 
than in the stagnation point itself, as proved from an 
analysis of the velocity distribution made by Homann. 



INTRODUCTION 



The prediction of the resistance of moving bodies 
is predicated on the knowledge of the conditions under 
which the laminar flow in the frictional layer adjacent 
to the wall becomes turbulent. The analytical treatments 
by Tollmien (references 1, 2), Schlichtung (references 3, 
4; 5, 6), and Gortler (references 7, 8) dealt with the 
two-dimensional problem. These. invest igat ions were based 
upon the method of small oscillations by ascertaining 
whether the small oscillations of various frequencies 
superposed on the basic flow are amplified or damped at 
given Reynolds numbers. A critical Reynolds number, the 
so-called "stability limit," can be considered as a first 
orientating measure of the stability of laminar flow, 
below which the oscillations of all frequencies for the 

velocity profile concerned are exactly damped.. This 

stability limit can be indicated for the velocity profiles 
accompanying pressure gradients and pressure rise in a 
two-dimensional frictional layer and the data of these 
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calculations used for an approximate prediction of the 
transitional inception on aircraft wings. 

An analysis of a surface of rotation in a flow along 
its longitudinal axi3, such as is approximately represented 
"by the airplane fuselage the resistance of which plays an 
increasingly important part in the power absorption of 
aircraft, would be desirable but for the following difficulty: 

The stability study requires the curves of the laminar 
velocity distributions to be known so accurately that even 
their curvatures can be reliably indicated; but at the 
present time there is little prospect of attaining this 
degree of accuracy for a surface of rotation of random 
sect ion. 

For this reason, the stability theory of the laminar 
frictional layer is to be explored on the simplest surface 
of rotation, the sphere, for which sufficiently exact ve- 
locity distributions in the frictional layer are likely 
to be available in the near future. 

In the concluding chapter;, the only reliable distri- 
bution known up to now, Homann's stagnation point profile 
(reference 9), is subjected to a stability analysis. 



II. GENERAL DIFFERENTIAL EQUATION OF THE DISTURBANCE 



The derivation of the stability theory for the sphere 
flow proceeds from the general Navi er-Stokes equations, 
expressed in spherical polar coordinates R, and $• 
(fig. 1) with u, v, and w denoting the velocities in 
direct ion, -of increasing 4 , R, and $, it affords for 
for the ; .axially. symmetrical problem 



The equation of continuity reads. 



1_ 9_ 

R sin * d 4 




(2) 



and the equations of motion for u and v have-, the form 
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The are length ■ starting 
from stagnation point and the 
distance n from the sphere 
surface are now introduced as 
coordinates 

s =R 0 », 
n = R — R t> 

with R Q as sphere radius. 

Restricted to the flow region 
wherein n «*; R 0 , the elimina- 
tion of pressure from (3) and 
(4) gives the equation of 
mot ion 



(4) 



lem then follows' the introduc- 
tion of the flow with Telocity 
components 

produced by superposition of 
basic flow and Interference 
motion, in equation (5). 

Restricted to the terms linear 
in cp and with those of subordi- 
nate order of magnitude discount- 
ed, the general differential 
equation reads for the present: 



a* u a_u — a a* u 

0(8n ' 5n5s dsdn 

u a v b*o a « a u 



8u a u . a*a . . v a_u 

"5Sb^ + %n' + R„ On 



R a !>n 

v d ' V 
bsdn 



dtds as 8s 



8'd 

*8s'" 



ft v 8 v 
8s 8n 



, 2 u8u f d* u . 8' a , 2 6' u 

^r, a* ~~ * U»« " 1 "a»ds«" t "fl 0 a n » 



, cotg »//?„ a* u i an 8'd 

" ' " dnds iJ 0 »sin«s/i? 0 8n dn'ds 



8*u cotga/i?,, 8'u 



8«« 

+ ii 0 a«« ^ 



« 0 ds» 1 fl„» sin' s/R„ d s 
2cotg s/i?, 8 a 
J?„» 8s 



8d , 2 8i) 
+ 8s 



'2k ] 
"iJo'sin's/ZtoJ 



■ (6) 



Let V and V designate the 
tangential and normal components 
of the axially symmetrical base 
flow about the sphere having 
stability to small vibrations is 
to be analysed. The superposed 
interference motion has the 
stream function 

v (,,n 1 t) = 9 .(s,r,) e 4a <'-' : ' ) . ...... (6) 

a indicating the spatial natu- 
ral frequency and the real part 
of c = c r + ici its phase ve- 
locity. The interference ve- 
locities are given by 



u* = - 



R, sin s/R, 8 n 

i a*> 



R, sin s/R a 
I 



,<«(«-<: I). 



' R, sin s/R, d s R, sin s/R, 



where 1 signifies differenti- 
ation with respect to n. 

The fundamental motion equa- 
tion of the interference pro b- 



(7) 



■V"<p+ -~{U<p'-U'<p) 



= - l i 1<p x v - 2 *« V " + «« ? i + i [ r (p'" - of ?.o - 

+ i[( 2 «. c _3a. t .-, l 7' 0 .|£+ P ^l] 
+ ^ iotg s/i} 0 [(3 V - c) a* y - 2 ( P ,>" + V »>')] (8) 

by referring the velocities to 
the velocity of undisturbed flow 
U^, the arc length to the sphere 
radius R Q and the wall distance 
as well as the wave length of the 
vibration to- a reference length 
of the friction layer, say, the 
di eplacement thickness 8* in a 
fixed point of the sphere surface 
equation (8) can also be written 
the form 

Expressed nondimensionally 

[U-c] {(p" — ot 2 <p) — U" <p + -^*- (U <p' — V'<p) 



l<p IV — 2<x 2 <p"-\--x.'<p\ 
-«*<p')-V"<p'] 



ot Reco 

+ f (2 «■ c - 3 «« V - V") ¥ + T^l 

ai?oL 8s 1 8sJ 



id* 

kR„ 



cotg s/R 0 I (3 V — c) n? <p — 2 ( U <p" + V <p') ] (9) 



This linear partial differen- 
tial equation of the fourth or- 
der represents the most general 
motion equation of the disturbed 
friction layer flow on the sphere. 
It differs from the corresponding 
equation of two-dimensional" flow 
by the addition of the underlined 
terms which have the same order 
of magnitude as the other term* 
on the right-hand side of equation 
(9) . 

The interference equation ap- 
plies equally in direct proximity 
to the stagnation point s « 0 
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because U and c decrease as 
sin s/R 0 , while cotg s/R 0 -* »; 



The normal component 7 of 
the velocity is obtained from (2), 



for R e 



.« it becomes, as expect- which after omission of the term 



ed, the two-dimensional interf er- §.3. simplifies to 



ence equation. 

The next problem is to prove 
that the interference equation 
for the flow past a sphere has 
exactly the same particular solu- 
tions as the two-dimensional, 
despite the additive underlined 
terms. 

III. SOLUTIONS 
Assuming a Re,* as very large 
and bearing in mind that 

£-■&.■'.■■•■■■•<«» 
the partial differential equa- 
tion (9) reduces to Tollmien's 
frictioniees interference equa- 
tion for two-dimensional flow 

(17 — c)(p" — a 1 ?) — V"(f = 0. . . . (11) 

in which "f rictionless solu- 
tions" qoj ,q> 8 are identical with 
those of the two-dimensional 
problem by reason of the same 
boundary conditions; hence all 
statements made elsewhere apply 
here also. In the construction 
of the solution cp e in the crit- 
ical layer U = c the friction 
must be taken into account in 
form of a transitional function, 
as explained in the following. 

The laminar basic velocity in 
the critical layer is approxi- 
mate to a parabola 

V-c=V,'{n-rn.)+¥f(n — n.)* . . . (12) 

with the new variables 

n = Ji= r^-' « = («**» tVT v '; f=«, • . (13) 



defined as in the two-dimension- 
al problem and 



1 dn, , r) dU.' 
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the interference equation (9) 
then reads: 
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(18) 



Limited to the terms linear in e 
because of the assumably large 
Reynolds number Re* equation 
(15) affords, after*multiplica- 
tion by e/u 0 ' 
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But then 

cC dtp , / d_p \ 

Of? "^iicco V <>V V l 

= Ccae*^+U.'«C*>[ V ^-^, (20) 

so that the underlined terms in 

(19) are disregarded. 

The friction correction for the 
frictionless flow cp 8 in the 
critical layer follows from the 

c-f- U.' eij 
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solution cp 2 in powers of e at 

9 j, = 9 t- -f- e-f u + «*»>»» + • • • .(.(HI 

which joins the frictionless so- 
lution ._ at some distance : from _ 
U = c. 

Posting (21) in (20) we get 
with cp 80 = 1 and multiplying 

.by -U 



1 di/ o"„, U„" dn e 

" r 6" ~dT 0 + U ° IT ~ "IT "ZT e 

- c c 4t _ i cotg 4/fl ° + T "•*) ] 

(22) 

• - . . 5u, bv 1 bp ' lb*v , 1 bo . b*v\ 

Again disregarding the terms with "^: + «^ = -jEJ^ , '\^ + Tb7 + bz'p 
€ the differential equation for 
the transitional function cp 8l 
is the same as in the two-dimen- 
sional pro hi em: namely. 



i* * = * — 



The Telocity distrihution in 
the stagnation point of the 
■sphere is to he analyzed as 
regards its stability against 
small vibrat ions. 

This velocity profile, which 
has the same form in the stagna- 
tion point of any other surface 
of ro.tation with other than van- 
ishing curvature radius, was 
derived by Hqmann (reference 9) 
from the hydrpdynamie equations 
of motion expressed for this 
purpose in cylindrical polar 
coordinat est 



• 0 u . 
dr bz 



b u 1 ftp j_ l b* u . 1 :5u 



0z> 



(25) 



br)< lV b n * - £/„' 
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Finally it can be readily proved 
that the frict ion- induced solu- 
tions in wall proximity cp 3 ,ep 4 
themselves are as in the two-di- 
mensional problem. For t -»» 0 
equation (19) affords the same 
differential equation 



where u the velocity along the 
increasing distance r from the 
stagnation point in the tangen- 
tial plane at the stagnation 
point and v the Telocity along 
the increasing distance % from 
the stagnation point past the 
body axis. 

Outside of the friction layer 
the potential velocity in the 
stagnation point is the same as 
that toward the plate 



(24) 



And since the limiting condi- 
tions of the differential inter- 
ference equation for flow past 
the sphere are also the same- as 
in the two-dimensional problem, 
it proves that the study of sta- 
bility of laminar flow past the 
sphere follows the very same 
method as in the two-dimensional 
p ro bl em . 

This result suggests the suk 
picion that the method of the 
two-dimensional problem is not 
merely applicable to the sphere 
but generally to any surface of 
rotation. Tor rigorous proof of 
this suspicion we would have to 
proceed from the equations of 
mo.tion in general orthogonal 
coordinates. 

It. THE STABILITY 07 LAUINAE FLOW 
IN THE STAGNATION POINT OF A 
SURFACE 'OF ROTATION 



u a " P r i v a = " 2 P * 

with t'he dimensionless wall 
distance 



(26) 



(27) 



and the formula; for the stream 
f unct ion 



= y p v r 8 f (t). 



(28) 



Homann's differential equation 
for the stagnation- point pro- 
file on the surface of rotation 
reads: 

f • + 2 f f« - f « 8 + 1 = 0 (29) 
with f« (J) » U - £il the dimen- 



sionless velocity in the fric- 
tion layer. 
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It differs from Hiemenz's differential equation for two- 
dimensional stagnation-point flow merely by the number 2 
in the second term; it has the same limiting conditions 
f = ft = o for t = 0 , and f« = 1 for £ — The 
solution is illustrated in figure 2. 

For the stability study it further requires an ac- 
ceptable approximation of the velocity distribution, for 
which in support of previous calculations (reference 10) 
an expression of the form 

E- = i _ (i ^ y ) n ( n = ,2,3,4 ...) (30) 
u a 

is employed, in conjunction with the equation illustrated 
in figure 2: 

JL = 1 - (1 - 0.4391 O 3 (31) 



Figure 3 contains the polar diagram for determining the 
neutral stability curve (fig. 4) which separates stable 
from unstable interference attitudes in the plane spanned 
by Reynolds number and interference frequency. The neu- 
tral stability curve for the two-dimensional stagnat ioa.- 
point profile according to Hiemenz is included for con- 
trast. It is seen that the stability limit Re *ritical 

on the three-dimensional stagnation-point profile is 
about a third of that of the two-dimensional profile. 
As a consequence, the laminar flow near the stagnation 
point of a surface of rotation is obviously more unstable 
in the case of small vibrations than the two-dimensional 
stagnation-point flow. 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Figs. 1,2,3,4 




Figure 1.- Definition of spherical polar coordinates, 
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Figure 2.- Velocity distribution in the friction layer at the stagnation 
point of a surface of rotation. 
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Figure 3.- Polar diagram for xu ■ LU " xu AU ~ • LU " A ° 

determining the neutral Figure 4.- Ueutral stability curve for 
stability curve. the stagnation point profile. 



